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ABSTRACT 

Mid-infrared photometry provides a robust technique for identifying active galax- 
ies. While the ultraviolet to mid-infrared (A ^ 5fim) continuum of normal galaxies 
is dominated by the composite stellar black body curve and peaks at approximately 
1.6 fiui, the ultraviolet to mid-infrared continuum of active galaxies is dominated by 
a power law. Consequently, with sufficient wavelength baseline, one can easily distin- 
guish AGN from stellar populations. Mirroring the tendency of AGN to be bluer than 
galaxies in the ultraviolet, where galaxies (and stars) sample the blue, rising portion 
of stellar spectra, AGN tend to be redder than galaxies in the mid-infrared, where 
galaxies sample the red, falling portion of the stellar spectra. We report on Spitzer 
Space Telescope mid-infrared colors, derived from the IRAC Shallow Survey, of nearly 
10,000 spectroscopically-identified sources from the AGN and Galaxy Evolution Sur- 
vey. Based on this spectroscopic sample, we find that simple mid-infrared color crite- 
ria provide remarkably robust separation of active galaxies from normal galaxies and 
Galactic stars, with over 80% completeness and less than 20% contamination. Con- 
sidering only broad-lined AGN, these mid-infrared color criteria identify over 90% of 
spectroscopically-identified quasars and Seyfert Is. Applying these color criteria to the 
full imaging data set, we discuss the implied surface density of AGN and find evidence 
for a large population of optically-obscured active galaxies. 
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1. Introduction 

The dominant sources of energy production in the universe are fusion in stars and gravitational 
accretion onto supermassive black holes. The tight correlation between nuclear black hole mass 
and bulge mass (e.g., Magorrian et al. 1998; Tremaine et al. 2002) implies the processes are 
intimately connected. However, identifying an unbiased census of black holes in the universe remains 
challenging, hampering our ability to fully probe this connection. For example, the spectral index 
of the X-ray background is significantly harder than the soft X-ray spectrum of bright, unobscured 
quasars, implying that there exists a large population of heavily-obscured AGN. Identifying the 
sources responsible for the (hard) X-ray background was one of the primary motivators for the 
Chandra X-ray Observatory, and Chandra has successfully identified a population of high-redshift, 
heavily-obscured, luminous type II quasars (e.g.. Stern et al. 2002a; Norman et al. 2002) as well as 
a significant population of apparently-normal galaxies with optical spectra dominated by starlight 
but hosting significant X-ray emission and implying the presence of an AGN (e.g., Barger et al. 
2001; Hornschemcicr et al. 2001; Stern et al. 2002b). 

We show here that the Spitzer Space Telescope (Werner et al. 2004) also provides a valu- 
able probe of AGN demographics. The obscuring dust which hides AGN from ultraviolet, optical, 
and soft X-ray surveys should be a strong, largely isotropic emitter in the mid- to far-infrared 
(A > 20fim). Furthermore, the different spectral energy distributions (SEDs) of stars and AGN 
produces different colors for normal and active galaxies in the mid-infrared (A ^ 10/um). In par- 
ticular, the composite black body spectra of the stellar population of normal galaxies produces 
an SED which peaks at approximately 1.6 fim, while quasars have a power-law SED, oc 
One of the early, traditional techniques to identify quasars relied on identifying sources with an 
ultraviolet excess, targeting sources whose SED does not plummet on the blue side of the stellar 
peak (e.g., Schmidt &; Green 1983). Though fruitful, this technique fails in two cases: (1) obscured 
quasars will preferentially suffer attenuation of their ultraviolet flux, and (2) high-redshift quasars 
disappear from the ultraviolet due to absorption from the Lya forest. 

By contrast, we seek to exploit the many- fold recent improvement in mid-infrared sensitivity 
afforded by Spitzer to target sources having SEDs that do not decline on the red side of the stellar 
peak. This provides an approach to identifying AGN that is relatively insensitive to extinction by 
dust or gas, and it has been demonstrated with a variety of instrumentation, including the Two 
Micron All-Sky Survey (Cutri et al. 2001; Glikman et al. 2004), the Infrared Space Observatory 
(Laurent et al. 2000; Haas et al. 2004), and the First Look Survey by Spitzer (Lacy et al. 2004). 

We present here a study of the mid-infrared colors of nearly 700 spectroscopically-confirmed 
quasars from the 9 deg^ Bootes field of the NOAO Deep Wide-Field Survey (NDWFS; Jannuzi & 
Dey 1999; Dey et al. 2005; Jannuzi et al. 2005). The mid-infrared colors are from the IRAC Shallow 

Survey (Eisenhardt et al. 2004), a guaranteed-time program with Spitzer. The spectroscopy comes 
from the AGN and Galaxy Evolution Survey (AGES; Kochanck et al., in preparation). AGES 
includes redshifts and spectral classifications for nearly ten thousand normal galaxies and nearly 
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one thousand AGN. Based on this spectroscopic sample, we show that mid-infraxed photometry 
provides a robust technique for identifying AGN. Using empirical color criteria to isolate AGN, we 
use the entire IRAC Shallow Survey to study the surface density of active galaxies, and evaluate 
the efficacy of mid-infrared selection of AGN. A coordinated paper discusses the mid-infrared colors 
of X-ray sources in this field (Gorjian et al. 2005). 

Section 2 briefly summarizes the two surveys used. Section 3 describes the mid-infrared prop- 
erties of the spectroscopically-identified sources, followed by a discussion of mid-infraxed selection 
of AGN in §4. The results are summarized in §5. Unless otherwise noted, all magnitudes refer to 
the Vega system. 

2. Survey Data 

2.1. IRAC Shallow Survey 

The Infrared Array Camera (IRAC; Fazio et al. 2004) is a four-channel instrument on the 
Spitzer Space Telescope which provides simultaneous broad-band images at 3.6, 4.5, 5.8 and 8.0 fim 
with unprecedented sensitivity. The IRAC shallow survey, a guaranteed-time observation program 
of the IRAC instrument team, covers 8.5 deg^ in the NDWFS Bootes field with three or more 
30 second exposures per position. Eisenhardt et al. (2004) presents an overview of the survey 
design, reduction, calibration, and initial results. The survey identifies ~ 270,000, 200,000, 27,000 
and 26,000 sources brighter than 5a limits of 12.3, 15.4, 76 and 76 //Jy at 3.6, 4.5, 5.8 and 8.0 
fim, respectively, where throughout IRAC magnitudes were measured in 6" diameter apertures and 
corrected to total magnitudes assuming sources are unresolved at the 1''66 — 1''98 resolution of 
IRAC (Fazio et al. 2004). The corresponding magnitude limits are 18.4, 17.7, 15.5, and 14.8 mag. 

2.2. AGN and Galaxy Evolution Survey (AGES) 

AGES (Kochanek et al., in preparation) is a wide-field redshift survey in the NDWFS Bootes 
field using Hcctospcc (Fabricant et al. 1998), a new, multi-object fiber spectrograph at the MMT 
Observatory. For the 2004 observing season, AGES targeted sources from the NDWFS catalog 
(Jannuzi et al. 2005) with R < 21.5, with priorities based on the multiwavelength photometry 
available. AGES observed (1) all extended sources with R < 19.2, (2) a randomly-selected 20% of 
extended sources with 19.2 < i? < 20, (3) all R<20 extended sources with IRAC 3.6, 4.5, 5.8 and 
8.0 /xm magnitudes < 15.2, 15.2, 14.7 and 13.2, respectively, and (4) sources with a strong selection 
bias in favor of counterparts to X-ray sources from Chandra X-ray Observatory imaging of the field 
(Brand ct al. 2005; Murray ct al. 2005; Kcntcr ct al. 2005), radio sources from the FIRST survey 
(Becker, White, & Helfand 1995), and 24 /xm sources with non-stellar J— [24] colors, where the 
24 /xm data derives from Spitzer Multiband Imaging Photometer (Rieke et al. 2004) observations 
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of the field (Le Floc'h et al., in preparation). Kochanek et al. (in preparation) provides a more 
detailed description of the target selection. In particular, point sources (in the Bw-, R-, or /-band 
images of NDWFS) were selected against unless they were associated with R < 21.5 counterparts 
to X-ray, radio, or 24 //m sources. 

We use the version 1.11 AGES catalog, containing all AGES spectra taken in 2004. The 
redshift catalog, created by N. Caldwell, P. Green, and C. Kochanek, includes 10,452 sources with 
confident redshifts. Non-stellar spectra were initially assigned to two rough classes: broad-lined 
AGN and galaxies. While all broad- lined sources correspond to obvious AGN, e.g., quasars and 
Seyfert Is (Syls), identified galaxies include both normal (i.e., apparently inactive) galaxies and 
narrow-lined AGN. Fairly conservative emission line diagnostics^ were then used to distinguish 
AGN from normal galaxies. The AGES sample with well-determined redshifts is comprised of 9394 
optically-inactive galaxies, 733 broad-lined AGN, 118 narrow-lined AGN, and 207 Galactic stars. 
Most of the stars are F stars, used for spectrophotometric calibration. Matching these sources 
with the 3.6 /xm-selected IRAC Shallow Survey catalog by requiring < 2" separation between the 
IRAC mid-infrared source and NDWFS optical source and > 5a detections in all four IRAC bands 
leads to a restricted spectroscopic sample of 4693 sources, comprised of 3959 normal galaxies, 576 
broad-lined AGN, 99 narrow-lined AGN, and 59 stars. For some analyses, wc only require IRAC 
5(7 detections in the two bluest IRAC bands; the matched catalog for this subsample contains 9432 
sources, comprised (as classified by AGES) of 8460 normal galaxies, 696 broad-lined AGN, 104 
narrow-lined AGN, and 172 stars. 



3. Mid-Infrared Properties of AGN 

Fig. 1 presents an IRAC color-color diagram for the AGES sample, with symbols denoting 

spectroscopic classification. Eisenhardt et al. (2004), based solely on photometric data, noted a 
vertical spur in the [3.6] — [4.5] versus [5.8] — [8.0] color-color diagram (see Fig. 5c of that paper) and 
suggested that it may be associated with AGN since a large fraction of those sources are spatially 
unresolved at 3.6 /xm. These spectroscopic results confirm that suggestion: mid-infrared colors 
clearly separate AGN from both stars and galaxies in the AGES sample. In order to achieve a 
physical understanding of Fig. 1, we now discuss the mid- infrared properties of stars, galaxies, and 
AGN. 

Galactic stars primarily reside in a restricted locus in mid-infrared color space, corresponding 
to the Rayleigh-Jeans tail of the blackbody spectrum: stellar photospheres have approximately 
zero color in the Vega system. The majority of the spectroscopically-confirmed stars in the AGES 



'^Objects fit with tlic galaxy templates wore flagged as AGN if they satisfled one of the following criteria. The 
first criterion was [NII]6583/Ha > 0.5, [OIII]5007/H/3 > 2, and detection of [OI]6300. The second criterion was 
the detection of [NeV]3426 emission. The first criterion can be used at lower redshifts {z < 0.35) while the second 
criterion can be used at higher redshift {z > 0.15). 
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sample confirm this expectation. The AGES sample also includes several stars which reside far from 

the stellar locus, showing mid- infrared excesses in the IRAC bands. Such stars likely host cooler 
circumstellar material in the form of post-asymptotic giant branch shells or circumstcUar disks. A 
more detailed treatment of the Galactic color-space outliers is deferred to a future publication. 

Because of its i?-band magnitude limit, AGES identifies normal galaxies out to z ~ 0.6. As 
seen in Figs. 1 and 2, galaxies at z < 1 span a relatively modest range in the bluest IRAC color 
combination, —0.2 < [3.6] — [4.5] < 0.5. For reference, we show the color evolution of two galaxy 
templates from Devriendt, Guiderdoni, k. Sadat (1999): M82 is a starburst galaxy with strong 
polycyclic aromatic hydrocarbon (PAH) features, while VGG1003 (NGC 4429), a Virgo cluster 
SO/Sa galaxy, is not actively forming stars. As the 3.3 jim. PAH feature shifts through the IRAC 
4.5 /im band at 0.2 ^ z < 0.5, galaxies appear redder in the [3.6] — [4.5] color combination. At 
z > 1, which lies beyond the redshift range in which AGES is sensitive to normal galaxies, the 
templates redden once again as the 1.6 jiva maximum in the photospheric emission from stellar 
populations shifts from the 3.6 jim. IRAC band to the 4.5 /xm IRAC band. 

Galaxies at z < 0.6 span a large range of color in the reddest IRAC color combination, 
< [5.8] — [8.0] < 3 (Figs. 1 and 3). These extremely red colors are caused by the 6.2 and 
7.7 jim. PAH features in actively star-forming galaxies shifting through the IRAC 8.0 /xm band. 
The 6.2 and 7.7 /im PAH features are much stronger than the 3.3 jim PAH feature, explaining 
why the observed range of [5.8] — [8.0] colors for normal galaxies is much more expansive than 
the observed [3.6] — [4.5] color range. The 3.3 /xm PAH feature later causes a modest, broad red 
bump in galaxy template [5.8] — [8.0] colors at z « 1.3. At 2; > 2, the 1.6 jum bump causes the 
galaxy templates to redden once again. As seen in Figs. 1 thru 3, the Devriendt et al. (1999) 
M82 and VCC1003 templates approximately bound the observed range of AGES galaxies in both 
IRAC color— color space and in IRAC color— redshift space. The narrow-lined AGN in the AGES 
sample appear in both the region of IRAC color— color space dominated by galaxies and the region 
dominated by broad-lined AGN, discussed next; the narrow-lined AGN sample apparently includes 
sources whose mid-infrared flux can be dominated by either stellar emission or emission associated 
with a powerful active nucleus. 

The isolation of broad-lined AGN in mid-infrared color-color space is the most dramatic feature 
of Fig. 1. While stars are generally restricted to a single locus of zero color and z < 0.6 galaxies 
define a horizontal swath in Fig. 1, broad-lined AGN create a vertical branch in the diagram. This 
characteristic property is easily understood. First, the lack of strong PAH emission in powerful AGN 
restricts their observed [5.8] — [8.0] colors. Second, since the A ^ 5 ;uni flux of AGN is dominated 
by power-law emission rather than a composite stellar spectrum which peaks at ~ 1.6 /im, the 
[3.6] — [4.5] color of AGN is significantly redder than that of low-redshift galaxies. The tendency of 
AGN to be redder than galaxies in the mid-infrared, where galaxies sample the red, falling side of 
the composite black-body stellar spectrum, mirrors the tendency of AGN to be bluer than galaxies 
(and stars) in the ultraviolet, where we are sampling the blue, rising side of the stellar spectrum. 



-6- 



The power-law nature of quasar SEDs is consistent with previous studies at optical wavelengths. 
Vanden Berk et al. (2001) presents a composite quasar spectrum derived from over 2200 Sloan 
Digital Sky Survey (SDSS) spectra. Fitting the continuum at rest-frame wavelengths AAq 1300 — 
5000 A with a power law, f^, oc Vanden Berk et al. (2001) finds an average spectral slope 
a = 0.44 lb 0.1, comparable to the a = 0.57 ± 0.33 derived by Pentericci et al. (2003) from a 
sample of 45 high-redshift SDSS quasars imaged in the near-infrared. IRAC observations of AGN 
from the AGES survey suggest a modest steepening of the average quasar spectrum at longer 
wavelengths. Considering the 696 AGES broad-lined AGN with robust (> 5a) data in the bluest 
IRAC passbands, the sigma-clipped average [3.6] — [4.5] color is 0.65 it 0.20 (Vega), corresponding 
to a spectra index a = 0.73 ± 0.84. Considering the 576 AGES broad-lined AGN with robust 
photometry in all four IRAC passbands, the sigma-clipped average [5.8] — [8.0] color is 1.06 ± 0.21 
(Vega), corresponding to a spectra index a = 1.13 it 0.59. We note that the reddening of AGN at 
z > 1 in observed [3.6] — [4.5] color and at z > 2 in observed [5.8] — [8.0] color is suggestive of the 
rest-frame 1.6 fim stellar bump contributing to the observed mid-infrared flux. Kochanek et al. 
(in preparation) presents a more detailed treatment of the composite AGES AGN spectrum in the 
mid-infrared, derived from a principal component analysis of the IRAC photometry. 



4. Mid-Infrared Selection of AGN 

We adopt the following empirical criteria to separate active galaxies from other sources in the 
AGES spectroscopic sample described above (shown as a dotted line in Fig. 1): 

([5.8] - [8.0]) > 0.6 A ([3.6] - [4.5]) > 0.2 • ([5.8] - [8.0]) + 0.18 

A ([3.6] - [4.5]) > 2.5 • ([5.8] - [8.0]) - 3.5, 

where A is the logical AND operator. The left boundary protects against fainter, higher-redshift 
galaxies, and is not critical for optically-bright samples such as AGES. The right boundary merely 
approximates the outer range of AGN colors. Considering Fig. 2, we caution that these color criteria 
may preferentially omit AGN at z ~ 0.8 and 2; ~ 2. Out of a total of 681 sources which reside in 
this portion of mid-infrared color-color space in Fig. 1, 522 (77%) are spectroscopically classified 
as broad-lined AGN, 40 (6%) are spectroscopically classified as narrow-lined AGN, 113 (17%) are 
spectroscopically classified as galaxies, and 6 (1%) are spectroscopically classified as stars. Out of all 
576 AGES sources spectroscopically classified as broad-lined AGN and having good 4-band IRAC 
data, 522 (91%) meet these mid-infrared color criteria. Of the 99 AGES sources spectroscopically 
classified as narrow-lined AGN and having good 4-band IRAC data, 40 (40%) meet these mid- 
infrared color criteria. Furthermore, less than 3% of the nearly 4000 AGES normal galaxies would 
be (mis-)classified as AGN using these mid-infrared criteria. We note that sources classified as 
normal galaxies by AGES based on their optical spectra may still host an active, supermassive black 
hole: Chandra has identified a large population of luminous X-ray sources (L2-iokeV ~ 10^^ ergs~^) 
with apparently normal optical spectra (e.g.. Stern et al. 2002b). Similar sources might be below 
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the existing X-ray limits for the Bootes field (Brand et al. 2005; Murray et al. 2005; Renter et al. 
2005), but still be correctly identified as AGN by Spitzer. 

Optical surveys find an increasing number of quasars as a function of flux — the SDSS finds 
15 quasars per deg2 to i* = 19.1 (Richards et al. 2002), the 2dF Quasar Redshift Survey finds 35 
quasars per dcg^ to bj = 20.85 (Groom ct al. 2004), and the COMBO-17 survey AGN luminosity 
function implies approximately 90 quasars per deg2 to = 21 (Wolf et al. 2003). Using the mid- 
infrared selection criteria on the full IRAG Shallow Survey photometric catalog (ver. 1.1) of 14,099 
non-saturated sources with > 5cr detections in all IRAG bands, there are 2014 AGN candidates over 
the 8.06 deg^ of the NDWFS with four-band IRAG coverage. This translates to an approximate 
surface density of 250 AGN candidates per dcg^ to a 8.0 fim. flux limit of 76 //Jy. For a spectral 
index a = 0.75 (§3), the mid-infrared depth of the IRAG images correspond to i? ~ 21. The 
significantly higher surface density of IRAG-sclcctcd AGN relative to optically-selected AGN of 
comparable mid-infrared flux implies a sizable population of obscured quasars with attenuation of 
their observed optical light. 

The relatively shallow {R < 21.5) limit of the AGES spectroscopic follow-up may bias our 
results: AGES only identifies galaxies up to z < 0.6, while the Devriendt et al. (1999) SEDs 
plotted in Fig. 1 encroach upon our AGN selection criteria at z f« 1.4. Some of the sources classified 
as AGN based on their IRAG colors are likely moderate-redshift galaxies, fainter than the AGES 
optical spectroscopic limits. However, our implied AGN surface density is below that predicted by 
models of mid-infrared AGN number counts. For Q,\ = 0.7, = 0.3, and Hq = 70kms~^ Mpc~^, 
Andreani, Spinoglio, &: Malkan (2003) predicts ~ 9 x 10^ AGN per deg^ with an 8 /xm flux higher 
than 100 /uJy, while Treister et al. (2004) predicts only ~ 1250 AGN per deg^ to the same flux 
limit. The large discrepancy is due to differences in the assumed evolution of the AGN luminosity 
functions. Since the IRAG Shallow Survey only detects ^ 3000 sources per deg^ with an 8 fim flux 
higher than 76 juJy, the Andreani et al. (2003) predictions are clearly high. 

5. Summary 

We analyze mid-infrared photometry from the IRAG Shallow Survey of 4693 sources with well- 
determined spectroscopic redshifts from version 1.11 of AGES. This sample includes 576 broad-lined 
AGN and 99 narrow-lined AGN. We find that mid-infrared photometry provides remarkably ro- 
bust separation of normal and active galaxies: simple mid-infrared color criteria identify over 90% 
of the spectroscopically-identified broad-lined AGN and 40% of the spectroscopically-identified 
narrow-lined AGN. Our selection criteria rely on a mid-infrared excess in AGN SEDs, mirroring 
the traditional ultraviolet excess technique which has long been exploited to identify quasars. Only 
17% of the AGES sources which fall within our mid-infrared color criteria for AGN are not spectro- 
scopically classified as AGN. Many of these sources likely indeed host a luminous, active nucleus, 
where either dust obscuration hides the signature of an AGN at observed optical wavelengths, 
or the AGN signature is drowned out by the stellar emission from the host galaxy (e.g., Moran, 
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Filippenko, & Chornock 2002). Simple mid-infrared selection criteria imply a surface density of 
250 AGN candidates per deg^ to an 8 //m flux limit of 76 //Jy, nearly triple the surface density 
inferred from optical surveys to this flux limit. AGN missing from optical surveys may be of the 
obscured, type II variety which are implied by unified models of AGN and are invoked by models 
of the X-ray background. We show that the Spitzer Space Telescope is a powerful tool for studying 
AGN demographics; Spitzer will help identify a less-biased sample of super massive black holes in 
the universe, allowing us to probe the interconnection of fusion-driven and accretion-driven energy 
production over cosmic history. 

Wc thank M. Ashby and J. Hora for carefully reading the manuscript, and DS thanks E. Trcis- 
ter for helpful discussion regarding models of AGN number counts in the mid-infrarcd. CSK thanks 
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made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, Califor- 
nia Institute of Technology, under NASA contract 1407. Support was provided by NASA through 
an award issued by JPL/Caltcch. Spectroscopic observations reported here were obtained at the 
MMT Observatory, a joint facility of the Smithsonian Institution and the University of Arizona. 
This work also made use of images and/or data products provided by the NDWFS, which is sup- 
ported by the National Optical Astronomy Observatory (NOAO). NOAO is operated by AURA, 
Inc., under a cooperative agreement with the National Science Foundation. 
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Fig. 1. — IRAC colors of spectroscopically-identified objects from the AGES survey of the Bootes 
field. Axes indicate both the Vega and AB magnitude systems. Spectral classification of sources 
is annotated in the upper left. The < z < 2 color tracks for two non-evolving galaxy templates 
from Devriendt et al. (1999) are illustrated; dark bull's eyes indicate z = 0. M82 is a starburst 
galaxy, while VCC1003 (NGC 4429) is an SO/Sa galaxy with a star formation rate approximately 
4000 times lower. The dotted line empirically separates active galaxies from Galactic stars and 
normal galaxies. 
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Fig. 2. — IRAC [3.6] — [4.5] color evolution of galaxies, AGN, and galaxy templates. Symbols 
are the same as in Fig. 1. The left-hand panel shows the full observed redshift range, highlighting 
the modest average color evolution of AGN. The right-hand panel shows z < 0.7, highlighting the 
color evolution of normal galaxies. Horizontal dotted lines illustrate the expected mid-IR color of 
power-law SEDs, f„ oc z/"", for 0.5 < a < 2.0. 
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Fig. 3. — IRAC [5.8] — [8.0] color evolution of galaxies, AGN, and galaxy templates. Symbols are 

the same as in Figs. 1 and 2. The left-hand panel shows the full observed redshift range, highlighting 
the modest average color evolution of AGN. The right-hand panel shows z < 0.7, highlighting the 
color evolution of normal galaxies. Horizontal dotted lines illustrate the expected mid-IR color of 
power-law SEDs, fi, oc 



